Instability of a Pseudo-Relativistic Model of Matter with 
Self-Generated Magnetic Field 

Marcel Griesemer*and Christian Tix 
Mathematik, Universitat Regensburg 
D-93040 Regensburg 



Abstract 

For a pseudo-relativistic model of matter, based on the no-pair Hamiltonian, we 
prove that the inclusion of the interaction with the self-generated magnetic field leads 
to instability for all positive values of the fine structure constant. This is true no 
matter whether this interaction is accounted for by the Breit potential, by an external 
magnetic field which is chosen to minimize the energy, or by the quantized radiation 
field. 



I Introduction 

The stability of matter problem concerns the question whether the minimal energy of a 
system of particles is bounded from below (stability of the first kind), and whether it is 
bounded from below by a constant times the number of particles (stability of the second 
kind). Stability of the second kind for non-relativistic quantum- mechanical electrons and 
nuclei was first proved in 1967 by Dyson and Lenard jl], g]. Since the new proofs of Lieb 
and Thirring, and Federbush in 1975 stability of matter is a subject of ongoing interest 
dealing with more and more realistic models of matter such as systems with a classical or 
quantized magnetic field included or with relativistic electrons (see || and the references 
therein). Stability with relativistic electrons is more subtle because of the uniform 1 /length 
scaling behavior of the energy, which holds for massless particles (high particle-energy limit). 
Then the minimal energy is either non- negative or equal — oo, so that stability of the second 
kind becomes equivalent to the statement that stability of the first kind holds for any given 
number of particles. We simply call this stability henceforth. 

This paper is about a pseudo-relativistic model of matter which is stable, but which 
becomes unstable when the electrons are allowed to interact with the self-generated magnetic 
field. The self-generated magnetic field may be described using either an effective potential 
(the Breit-potential), an external magnetic field over which the energy is minimized, or the 
quantized radiation field. In all these cases we find instability for all positive values of the 
fine-structure constant. In contrast to most other models, where the collapse of the system, 
if it occurs, is due to the attraction of electrons and nuclei Q ||, || [?J (there would be no 
collapse without this interaction), the instability here is due to the attraction of parallel 
currents. 

The model we study is based on a pseudo-relativistic Hamiltonian sometimes called 
no-pair or Brown-Ravenhall Hamiltonian describing N relativistic electrons and K fixed 
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nuclei interacting via Coulomb potentials. The electrons are vectors in the positive energy 
subspace of the free Dirac operator and their kinetic energy is described by this operator. 
For a physical justification of this model see the papers of Sucher S |), for applications of 
the model in computational atomic physics and quantum chemistry see Ishikawa and Koc 
flO[ , and Pyykko 113] . The Brown- Ravenhall Hamiltonian yields stability for sufficiently 
small values of the fine structure constant and the charge of the nuclei [li], [Tg|. ||] ; there 
are further rigorous results concerning the virial theorem |Tfl] and eigenvalue estimates fl8|| . 

We are interested in the minimal energy of this model when it is corrected to account for 
the interaction of the electrons with the self-generated magnetic field. This correction may 
be done for instance by introducing a external magnetic field V x A to which the electrons 
are then minimally coupled and whose field energy is added to the energy of the system. 
The field A is now considered part of the system and hence the energy is to be minimized 
w.r.t. A as well. The minimizing A for a given electronic state is the self-generated one 
(to avoid instability for trivial reasons the gauge of A has to be fixed). The energy of this 
system is unbounded from below if Na 3 / 2 is large, a being the fine structure constant, even 
if the vector potential is restricted to lie in a two parameter class {jAq(6x) : 7, S G R + } 
where Aq is fixed and obeys a weak condition requiring not much more than Aq ^ 0. This 
is our first main result. It extends a previous result of Lieb et al. || and is reminiscent of 
the fact that a static non-vanishing classical magnetic field in QED is not regular, in the 
sense that the dressed electron-positron emission and absorption operators do not realize a 
representation of the CAR on the Fock space of the free field Jl9| . 

Alternatively the energy-shift due to the self-generated magnetic field may approximately 
be taken into account by including the Breit potential in the energy. The resulting model 
is unstable as well. That is, the energy is unbounded from below if Na 3 / 2 is large, no 
matter how small a is. This is our second main result. It concerns a Hamiltonian that is 
closely related the Dirac- Coulomb-Breit or Dirac-Breit Hamiltonian, which is the bases for 
most calculations of relativistic effects in many electron atoms [^|, [l2| . We mention that for 
a = 1/137 the energy is bounded below if TV < 39 and unbounded below if N > 3.4 • 10 7 
(Theorem [| and Theorem |^) . 

A third way of accounting for the self-generated field is to couple the electrons to the 
quantized radiation field. From a simple argument using coherent states (Lemma |^) it 
follows that the instability of this model is rather worse than the instability of the model 
first discussed. 

As mentioned above the instability with external magnetic field was previously found by 
Lieb et al. ||. Our result extends their result and our proof is simpler. The model with 
Breit interaction corresponds to the classical system described by the Darwin Hamiltonian, 
which has been studied in the plasma physics literature (see |20| an the references therein). 
This classical model is thermodynamically unstable as well [}2l|| . 

In Sections || and fv we introduce the models with external magnetic field, with 
Breit potential, and with quantized radiation field and prove their instability (Theorem [l], 
Theorem ||, and Lemma ||). In Section III we also discuss dynamic nuclei for the model with 
Breit potential. There is an appendix where numerical values for stability bounds on Net 3 ' 2 
given in the main text are computed. 



II Instability with Classical Magnetic Field 

We begin with the model of matter with external magnetic field. For simplicity the electrons 
are assumed to be non-interacting and no nuclei are present. We could just as well treat a 
system of interacting electrons and static nuclei and would obtain essentially the same result 
(see Remark 4 below). 

Consider a system of iV non-interacting electrons in the external magnetic field V x A. 
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The energy of this system is 

N 

8 N (^, A) = ty, J] D^A)i>) + — J |V x A^)! 2 ^ 

where D fl (A) is the Dirac operator D(A) = a ■ (— iV + a 1 / 2 A(a;)) + /3m acting on the /x-th 
particle, and the vector ip, describing the state of the electrons, belongs to the Hilbcrt space 

N 

H N = f\ A+L 2 (R 3 ,C 4 ) 
A+ = X (^oo)(D(A = 0)), 



or rather the dense subspace 2)jv = Hat D i? 1 [(R 3 x {1, . . . ,4})^]. That is, an electron is 
by definition a vector in the positive energy subspace of the free Dirac operator. We will 
always assume the vector potential A belongs to the class A defined by the properties: 

i) V • A = 0, 

ii) A(x) — > as \x\ — > oo, 
Hi) / |V x A\ 2 < oo. 



Notice that 7ijv is not invariant under multiplication with smooth functions, in particular 
it is not invariant under gauge transformations of the states. It follows that the minimal 
energy for fixed A is gauge-dependent. It can actually be driven to — oo by a pure gauge 
transformation (see Remark 3 below) . To avoid this trivial instability we fixed the gauge of 
A by imposing conditions i) and ii). 

The constants a > and m > in the definition of D(A) are the fine structure constant 
and the mass of the electron respectively. In our units h = 1 = c, so that a = e 2 which is 
about 1/137 experimentally. We denote the Fourier transform of a function / by / or $(f) 
and use p or k for its argument rather then x or y. Our first result is 

Theorem 1. Suppose A 6 A is such that Re[e ■ A(p)] < in B(0,e) for some e £ R 3 and 
e > 0. Then there exist a constant Ca such that for all a > 0, m > and N > Ca&~ 3 ^ 2 

inf £ N (ip,jA(Sx)) = — oo. 

iPeD N ,\\rp\\ = l, 7,5eR+ 

Remarks. 

1. It is sufficient that A € A n L 1 and J R3 A(x)dx ^ 0, since A is then continuous and 
A{Q) 7^ 0. Thus we have instability for virtually all non-vanishing A E A. 

2. The smallness of Na 3 ^ 2 is not only necessary but also sufficient for stability (sec j| 
Section 4]). 

3. If the condition ii) that A vanishes at infinity (and thus the gauge fixing) is dropped 
there is instability even for N = 1 and the theorem becomes trivial. In fact for N = 1 
and A(x) = a =/= 0, £n=i(^,jA) — (ip,D(Q)i()) +ja 1 / 2 a J ip + (x)aip(x)dx which, as 
a function of 7, is unbounded from below for suitable ip € A+L 2 (R 3 ,C 4 ). 

4. The statement of the theorem also holds for the system of electrons and static nuclei 
with energy £n(i^, A) + a{ip, V c ip) where 

n k z n 1 K Z Z 
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if both N and Z K are bigger than Ca« 3 ^ 2 and if the energy is in addition min- 
imized with respect to the pairwise distinct nuclear positions R K . (see the proof of 
Theorem |). 

5. Quantizing the radiation field does not improve the stability of the system (see Sec- 
tion 



IV). 



The only way to restore stability we know is to replace Hn by the A-dependent Hilbert 
space 

N 

-Hn.a = A X(o,oo)OD(A))L 2 (K 3 ,C 4 ). 

Obviously £n{iP, A) > forjA € Hn.a- In fact even £jv(t/>, A) + a{tp, V c ip) is non-negative 
for Z K and a small enough (3) . 

Proof of Theorem [?]. We will only work with Slater determinants and the following repre- 
sentation of one-particle orbitals. If u e L 2 (R 3 ; C 2 ) then 

? ( E(p) + m \ l/2 ( u{p) \ 

with E(p) = \Jp 2 + m 2 , is the Fourier transform of a vector ip <G A + L 2 , and the map 
u i-> V, i 2 (K 3 ; C 2 ) -► A + L 2 (M 3 ; C 4 ) is unitary. 

It suffices to consider the case m = and find a Slater determinant ip = ip\ A . . . A tpisr 
and 7,(5 e M + such that £n(">P, r yA(5x)) < 0. In fact by the scaling ip i— > ^5, A 1— > 
defined by = 8~ 3 / 2 u fl (8~ 1 p) and A^a;) = SA(Sx) we can then drive the energy with 
m > to —00 because £(4>s, As, m) = 6£(ip, A, m/S) and £(ip, A, m/S) — > A, m = 0) 
for <5 — ► 00. 

Choice ofip. Let Q be the unit cube {p 6 M 3 |0 < Pi < 1}, u(p) = (xg(p),0) T , and 
e € R 3 an arbitrary unit vector. Set 

u M (p) =«(p- XN^e-nJ, /*=1,...,JV (3) 

where A is a positive constant to be chosen sufficiently large later on, and (n M ) M= i.. .;v C 1? 
are the AT lattice sites nearest to the origin, i.e., max /i=1 jv \ n y\ is minimal. We define 
ip = tpi A . . . A tp N by 



J_ / u„(p) \ w = _P 
s/2 \<r ■ UpU»(p)J ' P \p\ 



^(p) = ^(..:r^) ) (4) 



which is (^) for m = 0. Then ip £ 7ijv and (ip^, ip v ) — (u^, u v ) — <5^„. Notice that 

|p - \N 1/3 e\ < N 1/3 for all p e supp(u AI ) (5) 

at least for large N (see the appendix), i.e., in Fourier space all electrons are localized in a 
ball with radius N 1 / 3 and a distance from the origin which is large compared to the radius 
(since A will be large). 

Since ip = ipi A ... A ipN and m > we have 

JV n . 

£ N (iP ) A)=^2(iP fi) MiP fl ) + a 1 / 2 J2 / J^x)A{x)dx 

m=i m=i J (6) 

+ ^- [ |V x A(x)| 2 dcc 

87T 7 
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where J \i{x) = ip^ l (x)ai^ fl (x). By definition of 

J M (p) = i(27r)- 3 / 2 /" u*(fc - p) [er(u k ■ a) + (w fc _ p • £r)er] u^(k)dk. (7) 

Replace here by its defining expression and substitute (k — AA fl / 3 e — n M ) i— * fe. Since 

a; fc+AAr 1 / 3 e+n f , ~~ * e as A — * 00 an d since u has compact support, it follows that J^(p) 
converges to the current 

J (p) = e (2tt)- 2/3 y - p)u{k)dk (8) 

as A — > oo. More precisely | J M (p) — «/o(p)| < CA _1 | Jo{p)\ for A > Ao where Ao and C are 
independent of \x and AT. From Jo(p)|p| _1 , A(p)|p| € £ 2 it follows that 

ljl(p)A(p)dp = Jjl( P )A(p)dp + 0(X- 1 ), A^co. (9) 

After a scaling A i— > A<5 we may assume Re[e • A(p)] < in the support of Jo rather then 
in B(0,e), so that (g) is bounded from above by some —c\ < for A > Ao where c\ and Ao 
are independent of \x and N. Observing finally that 



(^, |V|^) = / \^M\Mdp < (A + (10) 



for all /x, we conclude 



S N (ip,jA) < (A + 1)A^ 4 / 3 - a 1 / 2 1 Nc 1 + 7 



2 c 2 



'1 



= (Ao + 1)A^ 4 / 3 - a^N 2 
4c 2 

which is negative for Na z / 2 large enough. At the end we inserted the optimal 7. □ 

The theorem has the obvious corollary 
Corollary 2. There is a constant C such that for all a > 0, m> and N > CoT^I 2 , 

inf £n(iP, A) = — 00. 

ip££) N ,\\i>\\ = l;AeA 

This result is due to Lieb, Siedentop, and Solovej ||. 
Remark. It is sufficient that C = 1.4 • 10 5 or that N > 3A ■ 10 7 for a -1 = 137, see the 
appendix. 

To conclude this section we compute min^ 6 _4 £^{ip, A). This will provide a link to the 
instability with Breit-potential discussed in the next section. To exhibit the A-dependence 
we write the energy as 

Sjfty, A) = e N {ijj, A = 0) + a l/2 J J{x)A{x) + — J |V x A(x)\ 2 dx, 

where J(x) is the probability current density associated with -0. Its functional dependence 
on tf> is not crucial here. A straight forward computation shows that the Euler-Lagrange 
equation for A is — AA = ina 1 / 2 Jt where Jt is the divergence free - or transversal - part 
of J. Comparing this equation with the Maxwell-equation for A in Coulomb gauge, which 
is DA — 47TO! 1 / 2 Jt, we find that the minimizing magnetic field is the self-generated one up 
to effects of retardation. Solving the Euler-Lagrange equation gives 

min A) = £ N (i>, A = 0) - £ / Jt \ X) Jt ^ V) dxdy. (11) 

aga 2 J \x-y\ 
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Ill Instability with Breit Potential 
III. A Static nuclei 

We now consider a system of TV (interacting) electrons in the external electric field of K static 
nuclei. There is no external magnetic field but a self-generated one which is approximately 
accounted for by the Breit potential. The energy is now 

N 

£ N {i ! ,R) = (ij,{Y J D li + a (V c -B))i>) (12) 



where 



N 



j , \ • ■ I \ , a H ' ( X fJ. X v) ® a v ' ( X fl x u) I , ,, 

B = l^^c — I + ~ — I ( 13 ) 

/J,<v 1 I- 



and V c is the Coulomb potential defined in (|l|) . H denotes the A-tuple (Hi, . . . ,Rk) of 
pairwise different nuclear positions and = D^{A = 0). As before ip belongs toDjy C Hn- 
The interaction — aB is usually derived from the corresponding interaction in the Darwin 
Hamiltonian by the quantization p/m at p2| or from QED: treating the interactions of 
the electrons with the quantized radiation field in second order perturbation theory leads to 
a shift of the bound state energy levels approximately given by —a(ip, Bip) p3| . Important 
for our purpose is that 



exchange terms J 2 J \x — y 

for any Slater determinant tj) = tp± A . . . A ipN of orthonormal functions ip^ (see the proof of 
Theorem g). 

We are interested in the lowest possible energy 

E n ,k = wi£ N (ip,R) 

where the infimum is taken over all if) € Dn with ||-0|| = 1 and all A"-tuples (Hi, . . . , Rk) 
with Rj 7^ Hfe for j ^ k. Our second main result is 

Theorem 3. There exists a constant C such that for all a > 0, m > 0, K G N and 
Zi, . . . , Zk G R+ 

En,k = -oo 

whenever N,^2Z K > Cmax(a -3 / 2 , 1). If ^2 Z^ > 1 it suffices that C = 5 • 10 4 or - when 
a' 1 = 137 - that N = Z K > 3.4 • 10 7 . 

Remarks. 

1. Similar as in Section H, V c and hence the condition on Z K may be dropped. Then 
there is instability for N > Cmax(a~ 3 / 2 , 1). It is for completeness of the model we 
keep V c in this section. 

2. Without B the energy is proven to be non-negative aZ K < 2/tt for all k and if a < 1/94 

fjl (see also ||). One expects however stability even for aZ K < 2 (- + ^) 1 a < 0.12 
[ 13| , |Tqj , which would cover the atomic numbers of all known elements. 



At least partly this theorem can be understood from Corollary |2j, Equation (|ll|) and 
Equation (ph. 
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Proof of Theorem [|. To begin with we prove (14). Let tp = tpx A. . . A^at with (tp^, ip u ) = 5^ 
and let J(x) — Ylu=i ^ ( x ) a i ; ^( x ) be the current density of ip. Note that Jr,i{p) — 



£j=i(*tf - ^)Jj(p) and that 

4tt p lPj I / 

With defined by 



^) = ^5>(% + ^h = ^(E' 



a ■ x 65 a • x 



it follows that 



~ / Jt \ X)Jt \ V) dxdy = \ ® ^> " 2/)^ ® Vv> 

•/ I 33 2/1 2 ^ 



(15) 



2 



which is equation ( |l4|) . Similar as in the proof of Theorem [j] it suffices to consider the case 
TO = and to find a Slater determinant ip = ip\ A . . . A tpN and nuclear positions R\ , . . . , 
such that £n{iPi R) < 0- 

Choice of the nuclear positions. A beautiful argument given in |3| show that, after moving 
some electrons or nuclei far away from all others 

V^)<e + ^l^J ^ dxdy 

for suitably chosen nuclear positions. Here e > is the (arbitrary small) contribution of 
the particles moved away. The second term can be dropped if Y," =l Z\ > 1. We use the 
inequality obtained in |14|j to estimate it from above and find 



1 N 

{ip,V c ip) <£ + const— ^(^,£>^). (16) 
n=i 

The number N of remaining electrons obeys N < ^2 Z K + 1 which is the reason for the 
assumption on ^2 Z K . Of course the choice of the nuclear positions depends on ip, which has 
not been specified yet. 

Define one-particle orbitals ip^ and currents J ^ and Jo exactly as in the proof of The- 
orem |l| with e being an arbitrary unit vector in R 3 . The convergence J \i{p) — > Jq(p) as 
A — > oo now implies that 



i r jtMJt(v) dxd N 2 

\x-y\ 

> c 2 N 



2 J \x~y~\ dxdV + ° (A ' 



(17) 



for A > Ao, where Ao and c 2 > are independent of N. 

To estimate the sum of exchange- and self-energy terms in (R~5J) notice that 



47T ^ 

{ippteipuiBix- y)V>i/®Vv) = I — \J^,t(p)\ \ d P> ( 18 ) 
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where J fJ/l ,(x) = tp*(x)ckip l/ (x). After writing J ^ u (p) as an integral in Fourier space in terms 
of and u v similar as in (Q) it is easily seen, using the support properties of and u v , 
that 

I J^t(p)\ 2 < \J^{p)\ 2 < H2ir)- 3 X (\p + n, L - n v \ < V3). (19) 

The TV balls B(n v , y/3), v = 1, . . . , N all lie in the ball B(0, N 1 / 3 ) and cover a given point 
at most, say, 4 3 = 64 times (replace the balls by cubes with side 2\/3). Therefore ( [l9| ) 
implies 

24 



£ \J»x(p)\ 2 < 192(27T)- 3 X (|p + n M | < iV 1 / 3 ) < — x (|p| < 27V 1 / 3 ) 

U— 1 

which in conjunction with (|lq ) gives 
1 

o £(Vv ® Vv, B(x - y)i> v ® Vv> < — ^ 4/3 - (20) 

Z » — ' 7T 

Rewriting the energy using ( |l5| ) and inserting the estimates ([l6|), (|Tc|), ( |l7| ) and ( |20| ) wc 
arrive at 

S N (ip,R) < ci(l + a)iV 4/3 -c 2 aiV 2 , c 2 > 

which is negative for TV > const max(a -3 / 2 , 1). This proves the theorem. □ 

For small N and small a there is stability. A similar result for the energy in Section 1 was 
proved in ||. 

Theorem 4. Suppose a < 1/94, max K Z K < IjnaT 1 andN-l < 2(2/-k+it /2){oT l -6r l ). 
Then En,k > 0. Inserting a — 1/94 and a = 1/137 we find stability for N < 39 and 
max Z K < 59. 

Froo/. Since B(x) < 2/|x| on C 4 <g) C 4 and l/|ac| < <5 _1 £) on A+i 2 (R 3 ;C 4 ) where 6 = 
2(2/n + it/2) p3 one has by the symmetry property of the states in TL^ 

Furthermore 

1 ^ 

^c>--V£» M on ^ (22) 

a * — ' 

for all a > with a max Z K < ^ and <5g < 1 /47 by || , where the number g of spin states 
may be set equal 2 H . Inserting d2l]) and (S3) in the energy proves the theorem. □ 



III.B Dynamic nuclei 

Making the nuclei dynamical would improve stability if their kinetic energy were the only 
term we added to (|l2]). However if the nuclei are relativistic spin 1/2 particles like the 
electrons and if the Breit-potential couples all pairs of particles, taking their charges into 
account, then the instability will actually become worse. 

Let us illustrate this for a system of N electrons and K identical nuclei of spin 1/2 and 
atomic number Z > 0. These nuclei are described by vectors in the positive energy subspace 



8 



of the free Dirac operator with the mass M > of the nuclei. To prove instability we adopt 
the strategy of the proof of Theorem ^ and thus assume M = and m = 0. As a trial- wave 
function we take 

if> = (ipi A ... A ip N ) <g> (0i A ... A 0k) 

where ^ is defined by equations (||) and ([|) and <p K is defined like -0k except that e and iV 
are replaced by — e and K respectively. It follows that in the limit A — > oo we get N + K 
(charge-) currents, the nuclear ones being larger than the electronic ones by a factor of Z 
but otherwise identical. The Breit interactions thus gives a negative contribution to the 
energy of order a(N + ZK) 2 . While the parallel currents of the N + K particles add up, 
the opposite charges of the electrons and nuclei cancel themselves. In fact for ip defined as 
above 



N 

(ii>, v c ip) <J2 dxdy 



z 



K 

2 E 

K,<<7 

K N 

K—l f-t— 1 

N(N- 1) 



\M*)\ 2 \MyW 

\x - y\ 



dR 1 dR 2 



dxdR 



l \MR2 



\Ri-R 2 \ 

\M*)\ 2 \MR)\ 2 



(23) 



72 K(K- 



\x-R 
1) 



NKZ 



(i + oix- 1 )) 



= [(KZ - N) 2 - KZ 2 - N] {1/2 + OiX' 1 )), 

where I is the limit of the above double integrals as A — > oo. Hence (ip, V c ip) is negative, 
e.g., if KZ — N and A is large. To achieve this in the static case we had to choose the 
nuclear positions properly. It is instructive to recall how this was done. The total energy 
is bounded from above by c^N 4 / 3 + if 4 / 3 ) - c 2 a(N + KZ) 2 , c 2 > 0, for N = KZ and A 
large, and is therefore negative for N = KZ large enough. 



IV Stability and Instability with Quantized Radiation 
Field 

Instability for the model with classical magnetic field implies instability for the model with 
quantized radiation field without UV-cutoff. In fact, for each classical magnetic field there 
is a coherent state of photons which reproduces the classical field as far as the energy is 
concerned. If an UV cutoff is introduced the relativistic scale invariance of the energy is 
broken and stability of the first kind is restored. The lower bound depends on the cutoff 
and goes to — oo as the cutoff is removed. 

The state of the system is now described by a vector ^ S Hn ® F where T denotes 
the bosonic Fock-space over _L 2 (M 3 ) ® C 2 , the factor C 2 accounting for the two possible 
polarizations of the transversal photons, and the total energy of ^ is 

N 

+ (1 <g> 

2 

/ dk\k\a[(k)a x (k), 



Hf = 



9 



where 



A(x) := / dk [e\{k)e lkx ® a x (k) + e x {k) e - lkx <g> a A (fc) 

A=l^ 

=: A + (aj) + A+(a;)* 



is the quantized vector potential in Coulomb gauge. The operators a\(k) and a\(k) are 
creation and annihilation operators acting on T and obeying the CCR 



[a\(kx), aUk 2 )] = SpvSfa - k 2 



[a»(fe 1 ),a»(fe 2 )]=Q 



where a A = a A or at, and the two polarization vectors e A (fe) are orthonormal and perpen- 
dicular to k for each k G M 3 . We use dk as a short hand for (27r) _3 ' 2 (2|fe|) _1 ' 2 dfc, and 
the subindex of a^, V M and fi^ indicates that these one particle operators act on the /^-th 
particle. While we used Gaussian units in Section || and HI we now work with Heaviside 
Lorenz units. 



Lemma 5. For each A c i G A PI L 
that 



there exists a vector 6 G T (coherent state) such 
£ q N ed (^<E>e)=£ N (^,A cl ) 



for allip G D 



jV- 



Proof. Pick A cl G A fl L 2 (M 3 ) and define ?7 A ( fe) = (|fc|/2) 1 / 2 e A (fc) • A cZ (fe) so that A d (x) 



Ati(x) + A+(x)* with 



Next set 



A ti(*) = E / dk Vx (k)e x (k) ( 
x=i J 



Hn) ■= i E / dfc k(fc)flA(fc) + Vx(k)a[(k) 



(24) 



A=l 



and = e in Wfl G JF. is called a coherent state, it is normalized and most importantly 
it is an eigenvector of all annihilation operators 



a x (k)e = Vx (k)e. 

,/From (|2~i|), ( |25| ) and the definition of r? A (fe) it follows that 



(25) 



and 



a^A+ix^i) ® 9 = (a M A+(x p ) <g> l) ^ ® 6 



(6,17/6)= | dfc|fc|5> A (fc)| 2 = i | dfcfc 2 |2 c; (fc)| : 



Inserting this in the energy proves the theorem. 



□ 



If an ultraviolet cutoff is introduced in the field operator A(x) then stability of the first 
kind is restored for all N and a certain range of values for a and Z K . This follows from | p3 , 
Lemma 1.5] and || Theorem 1]. 
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Appendix 

To obtain the numerical values for the constants in Corollary || and Theorem || we follow 
the proof of Theorem |^, up to a few modifications and explicitly evaluate the constants in 
this proof. 

The main modifications are that the two-spinor u is now defined in terms of the (nor- 
malized) characteristic function of the ball with radius 1/2 contained in the unit cube 
{p|0 < pi < 1} and that the 4-spinors ip2^—i are defined in terms of the i/^'s by in- 
terchanging the components of u, while n^^-i runs over the N/2 or - if N is odd - the 
(N + l)/2 lattice sites of 1? closest to the origin. The balls simplify the computation of 
Jo(p) and the double occupation n2 M -i = «2p reduces the kinetic energy. To begin with 
we note that the n unit cubes of the lattice Z 3 which are closest to the origin, all fit in a 
ball of radius 

" ,/3 (i)" 3 + ^ 

In particular the N/2 or (AT+l)/2 unit cubes containing the supports of the spinors ip^, fi = 
1, . . . , N all fie in the ball of radius bN 1/3 centered at AN 1/3 e where b = 1/2 if N > 1.2- 10 7 , 
b = 3/5 if N > 5 ■ 10 3 and b = V% if N > 1 (the ball of radius \/3n 1/3 contains never less 
than n lattice cubes). This replaces equation (jl]) and implies, together with equations (0) 
and (|), that 



|J» - Jo(p)\ < 



\Jo(p)\, 



X-b 

Using this and |J (p)| = l/2(2vr)" 3 / 2 (l - p) 2 (2 + p) one finds 



X > b. 



> N 2 
= N 2 



47T * — - 

dp—J (p)TJ (p) 



12b 
X-b 



Air --- 

dp— \J (p)\ 2 

p £ 



(A.l) 



1 



186 
X-b 
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35^' 



where T is the 3x3 matrix with the components = 5ij — PiPj/p 2 - This replaces (|l7|). 
We proceed to bound the self-energy and exchange terms. Inequality (fL9ft becomes 



\J^(P)\ 2 < 3(2^)- 3 X (|p + n M - n v \ < 1), 



(A.2) 



because the support of u now has diameter 1 not y/i. Since the N balls B( n^, 1) cover a 
given point at most 8 times (recall that now riiv-i = Ti2v) inequality (A.2) leads to the 
bound 



1 48 

i/f„,B(x - y)Vv ® Vv> < — WV 4/3 



(A.3) 
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improving (j20|). The kinetic energy is bounded by (A + b)N 4 / 3 and (ip, V c ip) < since 
^2 Z K = N and Z\ > 1. In conjunction with (|A.l[) and (|A.3|) this gives 



A + 6 + — ba - aTV 2 / 3 — 1 - - 

7T 707T V A — 



(A.4) 



For b = 1/2, a" 1 = 137 and the optimal A this is negative for N > 3.4 ■ 10 7 . For b = 3/5 
and a > arbitrary this is negative for N > C max(a -3 / 2 , 1) with C — 43859 where A was 
chosen to minimize C. This explains the numbers in Theorem ||. 

Now drop the term (48/7r)6a in equation (A.4) which was due to the exchange- and self- 
energy terms. By equation ( fll"| ) what we are left with is a upper bound for ini^^AeA £j^(ip, A). 
It is negative for b = \/3, the optimal A and a 3 ^ 2 N > 134'863, or for b = 1/2 the optimal A, 
or 1 = 137, and N > 3.4 • 10 7 . This explains the numbers in the remark after Corollary |2[ 
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